Phase transitions in two-dimensional model colloids in a 
one-dimensional external potential 

Florian Biirzle and Peter Nielaba 
Department of Physics, University of Konstanz, 18Jf.51 Konstanz, Germany 

(Dated: February 1, 2008) 

Abstract 

Two-dimensional melting transitions for model colloids in presence of a one-dimensional external 
periodic potential are investigated using Monte Carlo simulation and finite size scaling techniques. 
Here we explore a hard disk system with commensurability ratio p = V^as/{2d) = 2, where is 
the mean distance between the disks and d the period of the external potential. Three phases, the 
modulate liquid, the locked smectic and the locked floating solid are observed, in agreement with 
other experimental and analytical studies. Various statistical quantities like order parameters, their 
cumulants and response functions, are used to obtain a phase diagram for the transitions between 
these three phases. 

PACS numbers: 64.60.-1, 82.70.Dd, 05.10.Ln 
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Phase transitions in systems with reduced dimensionahty have been of interest at least 
since the days of Landau 1| and Peierls on whose ideas the work of Kosterhtz, Thouless, 
Halperin, Nelson, and Young (KTHNY) 

mm 

is based to a large extent. 
Since the formulation of their successful KTHNY theory, two-dimensional (2D) systems 
with periodic substates received more and more attention. In their groundbreaking work, 

n 

Chowdhury, Ackerson, and Clark [9], first investigated a 2D colloidal system under the 
influence of a one-dimensional (ID) periodic potential. This was achieved by interference 
of two laser beams, yielding a ID interference pattern. At a strong enough light intensity, 
crystallization of the colloidal suspension was observed, provided that the periodicity d of 
the periodic potential was chosen to be commensurate to the mean particle distance a^. 
For this phenomenon, Chowdhury et al. coined the name "laser induced freezing" (LIF). 
LIF is, qualitatively, due to the suppression of thermal fluctuations transverse to the ID 
periodic potential. Surprisingly, a remelting of the crystal at even higher light intensities 



was also observed in systems with short-range interactions between the colloids |10l. 11 11. Il2|. 
Consequentially, this reentrance scenario was named "laser induced melting" (LIM). This was 
explained as generic, caused by suppression of phonon fluctuations transverse to the potential 
troughs, leading to a decoupling of neighboring rows {l^. An analytical justiflcation for this 



explanation was subsequently given in 



13|. 



It should be emphasized that it was also possible to confirm this scenario by simulation, 
although some earlier simulations were inconclusive with respect to the reentrance phe- 
nomenon (for a more extensive discussion see [l^). In particular, Chakrabarti et al. [l^ 
recognized LIM even before its experimental discovery. Other Monte Carlo studies also veri- 
fied the reentrance behavior by considering hard and soft disks, as well as Lennard- Jones-like 
i.e.a.tio„ po.e„«als A,so^™l .e„o™aH.a.io„ ,oup studies on .Ms 



subject have been successfully performed 



20 



2l|. 



It turned out, however, that theoretical attempts based on Landau expansion and related 
mean field techniques failed to provide an explanation for the complete phenomenology (see 
[isl for a short review). The latter is believed to be due to the incorrect treatment of 
fiuctuations in such theories while it is known that, unlike in three dimensions, fiuctuations 
play a vital role in low dimensional systems and therefore cannot be neglected. 



This uncomfortable situation did not change until Radzihovsky, Frey, and Nelson 
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published their theory which is based on an elastic model of the triangular lattice. Within 



2 



this framework, all phase transitions considered so far can be explained by dislocation un- 
binding. As a further key result, the authors predicted also the existence of new unobserved 
phases. According to their analysis, the appearance of new phases and hence the complex- 
ity of the phase diagram depends, with given crystal orientation, only on the value of an 
integral quantity p which they called commensurability ratio. The latter is a generalization 
of the concept of commensurability and is formally defined by p = ^/3as/{2d), using the 
notation convention from For p = 1 this coincides with the previous meaning of com- 
mensurability. In this case, the theory predicts two distinct phases - the modulated liquid 
(ML), which mimics the geometry of the periodic potential, and the crystalline phase, which 
Radzihovsky et al. named "locked floating solid" (LFS). This name reflects the dual char- 
acter of this phase, since the colloids are unpinned along the potential minima, but pinned 
perpendicular to the minima. The appearance of those two phases and also the predicted 
shape of the melting curve are in good agreement with observations from most experiments 
and simulations discussed so far. 

The next difficult case is p = 2, of which we will present a Monte Carlo analysis in this 
work. For this case, Radzihovsky et al. predicted an additional phase, the "locked smectic" 
(LSm). This phase is characterized by a spontaneous symmetry breaking of the discrete 
translational symmetry present in the modulated liquid, with equal occupancy of each po- 
tential minima, while in the LSm only every pth minima is equivalently populated. The LSm 
exhibits, in contrast to the LSF, only short-range correlations between colloidal positions in 
adjacent troughs and therefore does not resist shear deformations for displacements along 
the potential minima. 

According to [l3|, both the transition from LFS to LSm and from LSm to ML should 



exhibit a reentrance behavior. Additionally, our investigation was stimulated by the work 



of Baumgartl et al. 



23| who recently observed the LSm phase at p = 2 experimentally. 



For our analysis, we consider a system of hard disks with diameter a. They interact via 
a pair potential defined by 

{oo, n,- < a 
(1) 
0, > a 

where is the distance between particles i and j. These particles are confined to a two- 
dimensional box with dimensions x Ly, with L^/Ly = ^/3/2. This system is subjected 
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FIG. 1: Two-dimensional pair correlation functions. In (a) an almost unperturbed liquid (Vq = 
0.01, p* = 0.88) is shown. At the potential strength Vq = 5 the other images show the (b) modulated 
liquid [p* = 0.83), (c) locked smectic [p* = 0.86), and (d) locked floating solid {p* = 0.90). 



to an external potential 

V{x,y) = Vosm(^2'^x^ (2) 

which is periodic in the x-direction and constant in the y-direction (note that in [l^ the 
periodicity is chosen to be in the y-direction). The constant do in eq. ([2]) is defined as 
do = a/3/2 to meet the requirements from commensurability. Our system is completely 
characterized by two quantities, namely by the reduced density g* = ga"^ and the reduced 
potential strength V^* = Vo/iksT)^ where ks is the Boltzmann constant and T the temper- 
ature. For simplification, a was set to unity in our simulation. 

In this work we use standard Monte Carlo techniques In particular, the Metropolis 



algorithm 



25| is applied to the canonical (NVT) ensemble. Further we use the technique of 



block analysis, where we divide our system in several smaller subsystems and calculate the 
quantities of interest therein. This method allows us to compute many different system sizes 
within one simulation run. In this paper, we distinguish these subsystems by attaching an 
index L on the appropriate quantities. The value of L is calculated via S^/d with being 
the side length in the x-direction of a subbox. 
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All simulations considered here have been carried out with = 1024 particles, in order 
to compare the results with {l^. We executed our calculations by setting up an ordered 
crystal and reducing g* at fixed V^*. For equilibration, especially at high potential strengths, 
one must ensure that the system under consideration cannot be trapped in local free energy 
minima. In experiments, this problem can be solved by switching the laser off and on while 
in simulation this is of course not possible. To overcome this issue, we use nonlocal "trough 



moves", already introduced in 



1J|, in addition to the ordinary Monte Carlo moves. In 



these trough moves, that are used in every simulation run, it is attempted to place particles 
in other potential troughs. It was possible to show that such moves are required to 
reach equilibrium, since at large Vq the formation of dislocations is artificially hindered by 
considering only local moves, because the particles cannot bypass each other. 

Our simulations were carried out mainly on contemporary personal computers with 4x 10^ 
to 10 X 10^ Monte Carlo steps (MCS) from which 2 x 10^ to 4 x 10^ were used for relaxation, 
depending on statistical inefficiency. Typical simulation times ranged from about 50 to 120 
CPU hours. 

In order to distinguish all appearing phases in our simulation by visual inspection, we 
computed two-dimensional pair correlation functions. For systems with densities below crys- 
tallization and without external potential, the graphs of these functions should exhibit the 
typical concentric circles of an unperturbed liquid. With an external potential, one expects 
to see lines along the potential minima characterizing the modulated liquid. However, in 
the crystalline (LFS) case, the typical structure of the triangular lattice should crop up. In 
the LSm phase, compared to the modulated liquid, every second line must vanish due to 
breaking of translational symmetry. 

A survey of such 2D correlation functions at different V^* and g* already revealed some 
interesting properties of the phases. For potential strengths Vq* ^ 1, we found for densities 
below spontaneous solidification no significant disturbance by the external potential and 
hence no smectic phase, as one would expect for almost unperturbed liquids [see Fig. [I](a) 
as typical example]. At intermediate potential strengths with Vq > 1, we found a completely 
different behavior in the liquid regime. To exemplify this, we consider a potential strength 
Vq" = 5 at various densities g*. At g* = 0.83 we recognize a modulated liquid where every 
potential minima, on the average, is occupied with particles [Fig. [11(b)]. For a slightly higher 
density at g* = 0.86 we find indeed the predicted LSm phase, where only every second lane 
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is occupied [Fig. [11(c)]. Finally, at g* = 0.90 the structure of the triangular lattice can be 
observed, indicating the presence of the crystalline LFS phase [Fig. [11(d)]. We would like to 
mention, that these graphs are qualitatively in good agreement with Baumgartl et al. 
who obtained similar results for the appropriate phases. For larger potential strengths, it is 
recognized that the density range in which the LSm phase appears becomes both narrower 
and shifted to larger densities. We interpret this observation as a strong evidence for a 
remelting transition in this regime. 

Since the circumstances in the LFS phase are very similar to the appropriate phase in 
the case p = 1, we concentrate now on the phase transitions between ML and the new LSm 
phase. These freezing and remelting processes can qualitatively be described as follows. 
At low potential strength (V^* <^ 1) the ML is quite disordered, i. e. the particles are 
randomly distributed between the potential minima. The increase of the potential strength 
{Vq > 1) then leads to a reduction of fluctuations perpendicular to the potential troughs. 
Subsequently the particles occupy every second potential minima and form a quasi-long- 
range orientational order along the troughs. This is in contrast to the LFS, where also 
a quasi-long-range positional order exists; but it must be noted that still some particles 
or even particle groups may overcome the potential barrier and occupy free positions in 
adjacent troughs. This prevents complete decorrelation of neighboring particle rows; but 
at even higher potential strengths (V^q* ^ 1) the fluctuations perpendicular to the troughs 
are further reduced. Consequentially, particles can occupy adjacent rows without geometric 
restrictions due to particles from other rows. From this it can be concluded, that the entropy 
of the system in the ML phase would be higher than in the LSm phase. Since interaction 
energies are for large potential strengths approximately the same in both phases, the free 
energy in the ML phase becomes lowered and so a transition from LSm to ML is possible. 

For numerical calculation of the phase transitions points, we introduce two different order 
parameters. Both are Fourier transforms of the particle densities in direction of reciprocal 
lattice vectors. We consider only the set of the six smallest reciprocal lattice vectors of the 2D 
triangular lattice. To detect the locked floating solid phase, the vector Gi = 27r/(io(l/2e^. + 
\/3/2ey), enclosing an angle of 7r/3 with the wave vector K = An/doe^ of the external 
potential, is used. With Gi the definition of the appropriate order parameter t/'Gi reads 

N 

tpGi = J]^exp(-iGirfc) , (3) 

k=l 
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where is the position of particle k. Note that this definition was first introduced in 
\l\ for the treatment of the case p = 1. For the transition to LSm, we choose the vector 
Go = 27i/doex which is parallel to K and has half of its magnitude. The corresponding order 
parameter tpGo is then defined analogous to ^/'Gi- 

From these order parameters, the phase transition points have been calculated with the 



cumulant intersection method 



27|. The fourth order cumulant Ul is defined via 



uUq\K) = i-^j^ (4) 

where the index i has either the value or 1, depending on the appropriate order parameter. 
Since in continuous phase transitions the correlation length ^ diverges, the cumulant Ul = 
Ui^Las/i) becomes independent of the system size. From this it follows that cumulants for 
different system sizes intersect at the transition point. Note that there is also an intersection 



point in first order phase transitions [28| which makes it unnecessary for us to judge whether 
the transitions observed are of first order or continuous. 

It must be noted that, at high V^*, the cumulants to both order parameters do not 
show a true intersection point any more. Instead, there is a unification point, after which 
all cumulants collapse onto a single curve. This behavior is known to be typical to the 

n n 

anisotropic XF-model [29|, which was also observed in the study by Strepp et al. In 
those cases, at this unification point the phase transition was assumed to take place. 
Additionally, to supplement our cumulant analysis, we consider the response functions 

fc^TxG. = K^^j^ - (^G,)y , (5) 

also defined as in jl^. These susceptibilities are known to increase with increasing L, 
yielding a maximum in the density range where the phase transition occurs. Those maxima 
are shifted to lower densities compared to the cumulant intersection points, due to finite size 
effects {l^. Nevertheless, the form of the melting curve should be in agreement with those 
obtained with the cumulant intersection method. In our calculations we observe that this is 
indeed the case. For the sake of clarity, since a comparison with the configurations and 2D 
pair correlation functions shows that the results obtained by cumulant intersection points 
are more accurate, the susceptibility maxima have been omitted from the phase diagram. 

The phase diagram shown in Fig. [2] was obtained by using the data from the cumulant 
intersection points. As the most important result, we see that the melting curves for both 
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FIG. 2: Phase diagram in the p* /Vq plane. Transitions points have been obtained by considering 
cumulant intersection points. Open circles: order parameter ^Gi, and closed circles: V'Go- 



order parameters show a distinct remelting behavior at higher Vq, as was expected by the 
theory of Radzihovsky et al. By considering the curve composed of cumulant intersection 
points belonging to T/^n , (open circles), we find that the melting curve resembles those 



obtained by Strepp 



14| for p = 1 quite well. The other transition curve from ML to LSm 



(closed circles) shows that here the global minimum of the curve is slightly shifted to higher 
potential strengths. Also the minimum is located at considerable lower densities. Finally, 
it must be emphasized that at Vq -^0, the different melting curves collapse into one single 
curve, as is expected for physical reasons. 

In conclusion, we have investigated phase transitions in 2D model colloids by Monte Carlo 
simulation. In particular, melting transitions in the presence of a ID periodic potential were 
studied for the commensurability ratio p = 2. In contrast to the case p = 1, an additional 
intermediate phase between the LFS and the ML, the LSm, has been confirmed. This is 
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qualitatively in good agreement with theoretical predictions and experimental results alike. 
Furthermore, by defining appropriate order parameters, we were able to construct a phase 
diagram for the three phases observed. This is a substantial progress over the experimental 



work 



23| . since in that work only isolated points of the whole parameter space were captured. 



Finally, we would also like to encourage further experimental work on this exciting sub- 
ject, especially the construction of a phase diagram for this case experimentally. In addition, 
explorations of cases with even higher commensurability ratios, where further phases should 
appear, would be of great interest. 

We want to thank our former colleague Wolfram Strepp for his preliminary work on the 
subject and for providing us with the necessary details. This work was, in part, supported 
by the Deutsche Forschungsgemeinschaft (SFB TR6/C4). Granting of computer time from 
HLRS, NIC, and SSC is gratefully acknowledged. 
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